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ABSTRACT 



Context. Calcium is a key element for constraining the models of chemical enrichment of the Galaxy. 

Aims. Extremely metal-poor stars contain the fossil records of the chemical composition of the early Galaxy and it is important to compare Ca 
abundance with abundances of other light elements, that are supposed to be synthesized in the same stellar evolution phases. 
Methods. The NLTE profiles of the calcium lines were computed in a sample of 53 extremely metal-poor stars with a modified version of the 
Q ' program MULTI, which allows a very good description of the radiation field. 

5— I , Results. With our new model atom we are able to reconcile the abundance of Ca deduced from the Cal and Call lines in Procyon. This 
C/3 1 abundance is found to be solar. 

' -We find that [Ca/Fe] = 0.50 ± 0.09 in the early Galaxy, a value slightly higher than the previous LTE estimations. 
-The scatter of the ratios [X/Ca] is generally smaller than the scatter of the ratio [X/Mg] where X is a "light metal" (O, Na, Mg, Al, S, and K) 
I ' with the exception of Al. These scatters cannot be explained by error of measurements, except for oxygen. Surprisingly, the scatter of [X/Fe] 
^ , is always equal to, or even smaller than, the scatter around the mean value of [X/Ca]. 
^\ . -We note that at low metallicity, the wavelength of the Cal resonance line is shifted relative to the (weaker) subordinate lines, a signature of 
the effect of convection. 

t-H . -The Ca abundance deduced from the Cal resonance line (422.7 nm) is found to be systematically smaller at very low metallicity, than 
• the abundance deduced from the subordinate lines. Our computations of the effects of convection (3D effects) are not able to explain this 
■^-j- difference. A fully consistent 3D NLTE model atmosphere and line formation scheme would be necessary to fully capture the physics of the 
f*^) , stellar atmosphere. 

(N ■ 

■ Key words. Line : Formation - Line : Profiles - Stars: Abundances - Stars: Supernovae - Galaxy evolution 
•i-H ' 

X: 

£3 | 1. Introduction mass stars have been formed at the very early phases of the 

Galaxy and the chemical composition of their atmosphere re- 
An homogeneous sample of 53 metal-poor stars, most of them flectfj the yieMs of the first massive type n supe rnovae which 
extremely metal-poor (EMP stars with [Fe/H] < -2.9), has have a very short i ife _ time . These supernovae produce more 
been observed by Cayrel et al. d20Ql, and Bonifacio et al. << a _ e l e ments" (O, Mg, Si, S, Ca) than "iron-peak" elements. In 
42007] H00i. The aim of this paper is to determine more pre- contrast? less maS sive type I supernovae, which have a much 
cisely the calcium abundance in these stars. 120-125 F6.2 — longer life . time and explode laterj produce more iron .p e ak el- 
Coit18 NLTE corrections, model (7000,2,-3.5,0.3) These low ements than a . e i e ments. As a consequence in the atmosphere 
of the EMP stars formed at the beginning of the Galaxy a rel- 



Send offprint requests to: M. Spite ative enhancement of the a-elements (compared to the Sun) is 

e-mail: Monique.Spite@obspm.fr observed. The level of this overabundance is one of the fun- 

* Based on observations obtained with the ESO Very Large damental parameters of the chemical evolution models of the 

Telescope at Paranal Observatory, Chile (Large Programme "First Galaxv 
Stars", ID 165.N-0276(A); P.I.: R. Cayrel). 

** The NLTE collections of the Ca lines are available in elec- On the other hand, the relative production of the a-elements 

tronic form at the CDS via anonymous ftp to cdsarc.u-strasbg.fr in a supernova d epends on the mass of this supernova (e.g. 

(130.79.128.5) or via http://cdsweb.u-strasbg.fr/cgi-bin/qcat7J/A-i-A/ Kobayashi et al., 2006): [Ca/Fe] tends to be lower for more 
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massive supernovae. Therefore the level of [Ca/Fe] in the early 
Galaxy constrains the IMF in the early times. 

Moreover, Cayrel et al. (120041 and Bonifacio et al. (120091 
have shown that for [Fe/H] < -2.7 the abundances of the 
a-elements relative to iron ([Mg/Fe], [Si/Fe], [Ca/Fe], and 
[Ti/Fe]) are constant with a small scatter but that surprisingly 
the scatters of these elements relative to magnesium (another 
a-element) are larger. This larger scatter cannot be explained 
by a lower precision of the abundance of magnesium. Cayrel 
et al. ( 120041) and Bonifacio et al. ( 120091 ) computed the abun- 
dances under the LTE hypothesis. The non-LTE (NLTE) abun- 
dances of two a-elements Mg and S have been then computed 
in Andrievsky et al. (2010) and Spite et al. (1201 It . These papers 
confirm that the scatter of abundance ratios is generally larger 
when Mg replaces Fe as a reference element. 120-125 F6.2 — 
Corrl8 NLTE corrections, model (7000,2,-3.5,0.3) It is then in- 
teresting to check wether the link between the abundance of Ca 
and the abundance of the other a-elements is closer. 

In this paper we have carried out a NLTE analysis of the 
calcium abundance in the atmosphere of these EMP stars and 
have tried to estimate the influence of convection (3D compu- 
tations). 

For this new analysis we have used the resonance line of 
Cal and about 15 subordinate lines. We have also used the line 
of the Ca II infrared triplet at 866.21nm. Unfortunately, the two 
other lines of this triplet (at 849.81 and 854.21nm) are outside 
the observed spectral range. 

2. Star sample and model parameters 

The spectra of the stars investigated here have been presented 
in detail in Cayrel et al. (2004) and Bonifacio et al. (|2007T >. The 
observations were performed with the high-resolution spectro- 
graph UVES at ESO-VLT (Dekker et al.. l2000l . The resolving 
power of the spectrograph is R ss 45000, with about five pixels 
per resolution element and the S/N ratio per pixel is typically 
about 150. 

The fundamental parameters of the models (effective tem- 
perature T e ff, logarithm of the gravity logg, and metallicity) 
have been derived by Cayrel et al. (120041 ) for the giants and 
Bonifacio et al. ( 120071) for the turnoff stars. Briefly, tempera- 
tures of the giants are deduced from the colors with the cal- 
ibration of Alonso et al. (1999 2001), and temperatures of 
the tumoff stars from the wings of the Ha line. Moreover we 
checked that these Ha temperatures agreed with the tempera- 
tures derived from the color V-K and the calibration of Alonso 
et al. (1996). The gravities are derived from the ionization equi- 
librium of iron (under the LTE approximation) and we note that 
they might be affected by NLTE effects. The parameters of the 
models are repeated in Table Q~]for the reader's convenience. 

3. Determination of the calcium abundance 

The NLTE profiles of the calcium lines were computed with a 
modified version of the code MULTI (Carlsson, 1986 Korotin 
et al., 119991 ), which allows a very good description of the ra- 
diation field. This version includes opacities from ATLAS9 



(Kurucz, 1992), which modify the intensity distribution in the 
UV. 

3.1. Atmospheric models 

For these computations we used Kurucz models without over- 
shooting (Castelli et al., 1997). These models have been shown 
to provide LTE abundances very similar (within 0.05 dex) to 
those of the MARCS models used by Cayrel et al. (2004) 
and Bonifacio et al. ((2)09). The solar model was taken from 
Castellfl with a chromospheric contribution from the VAL-3C 
model of Vernazza et al. ( 119811 ) and the corresponding micro- 
turbulence distribution. 

3.2. Atomic model 

Our model atom of calcium is similar to the one used by 
Mashonkina et al. (2007) but it includes some more levels and 
more recent atomic data. Seventy levels of Ca I, thirty-eight lev- 
els of Ca II, and the ground state of Ca III were taken into ac- 
count; in addition, more than 300 levels of Ca I and Ca II were 
included to keep the condition of the particle number conser- 
vation in LTE. The fine structure was taken into account for the 
levels 3d 2 D and 4p 2 P° of Call. 

The energy levels were taken from the NIST atomic spec- 
tra database (Sugar & Corliss, 1985 ). The ionization cross- 
sections were taken from TOPBASE. The oscillator strengths 
of the Cal and Call lines were taken from the most re- 
cent estimations: Wiese et al. (119691 ). Smith & Raggett (119811 ). 
Smith (fT988l . Theodosiou (fT989l . Morton JT99TT ). Kurucz 
(1993), and from TOPBASE for the lines occuring between 
non-splitted levels. For the forbidden transitions, we used 
TOPBASE and Hirata & Horaguchi JT995I 

We considered 35 1 transitions in detail; for 375 weak tran- 
sitions the radiative rates were fixed. Collisional rates between 
the ground level and the ten lower levels of Ca I were taken 
from Samson & Berrington (2001). For Call, collisional rates 
were taken from Melendez et al. (12007) instead of those of 
Burgess et al. (1995) used by Mashonkina et al. (2007) for the 
lower seven terms. 

For the other transitions (without data) we used for allowed 
transitions, the Van Regemorter ( 1962 ) formula, and for the for- 
bidden transitions the Allen ( 1973 ) formula. 

Electron impact ionization cross-sections were calculated 
by applying the formula of Seaton dl962| ), with threshold pho- 
toionization cross-sections from the Opacity-Project data. 

Collisions with hydrogen atoms were computed using the 
Steenbock & Holweger ( 1984) formula. The cross sections cal- 
culated with this formula were multiplied by a scaling factor 
Sh- the "efficiency" of the hydrogenic collisions. This factor 
was constrained empirically by comparing the Ca abundance 
obtained from different lines in the Sun and in some reference 
stars (Procyon, HD 140283, HD 122563). The best agreement 
was obtained with Sh =0.1 which agrees well with Ivanova et 
al. (12002] ) and Mashonkina et al. d200"7j ). 



1 http://wwwuser.oats.inaf.it/castelli/sun/ap00t5777g44377klasp.dat 
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Table 1. Program stars and their parameters. The stars marked with an asterisk are carbon-rich. Columns 2 to 5 give the main 
parameters of the stars. Column 6 is the NLTE calcium abundance deduced from the subordinate Ca I lines with the number of 
lines and the standard deviation (columns 7 and 8). Columns 9 and 10 list the abundance of Ca derived from a NLTE computation 
of the 422.67 and 866.22nm lines. The last two columns give [Ca/H] and [Ca/Fe] based on the calcium abundance deduced from 
the Ca I subordinate lines. 
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Fig. 1. Profiles of Ca I and Ca II lines in Procyon. The wave- 
lengths are in A. The small crosses represent the observed spec- 
trum and the (red) line the computed profile. With our model 
atom, the observed and computed profiles of the Ca I and 
Ca II lines agree well. 



3.3. Consistency check 

To test the Ca atom model, we computed the profiles of the 
calcium lines in the Sun, in Procyon, and in two classical 
metal-poor stars HD 122563, and HD 140283. The synthetic 
spectra were computed following the procedure described in 
Korotin (2008): we calculated the departure coefficients factors 
"b" for the Ca lines with MULTI and then used these factors in 
the LTE synthetic spectrum code. 

• For the Sun we used the solar atmosphere model com- 
puted by Castelli with a chromospheric contribution following 
Vernazza et al. (I198U (see section 13. U . A micro turbulence 
velocity £ t = 1.0 kms -1 was adopted in the atmosphere. 

We computed the profile of 49 lines of Cal and 17 lines 
of Ca II. The log gf value of the lines and the broadening pa- 
rameter due to collisions with hydrogen atoms, logyvw/NH 
(for T=10000K), are given in Table We found loge(Ca) = 
6.31 + 0.05 from the Cal lines and loge(Ca) = 6.30 + 0.07 
from the Call lines. These values agree well with the mete- 
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Fig. 2. Profiles of Ca I and Ca II lines in HD 140283. The wave- 
lengths are in A. The symbols are the same as in Fig. Q] The 
observed spectrum and the synthetic profiles computed with 
log e(Ca) -4.12 agree well. 

oritic calcium abundance loge(Ca) = 6.31 + 0.02 (Lodders et 
al.,[2009) and the photospheric abundances derived by Asplund 
et al. ( 120091 1 for the solar atmosphere: log e(Ca) = 6.34 + 0.04 
(a value that includes 3D effects). 

For most of the lines in Table [2] the broadening parameter 
logyvw/Nn (for T=10000K) due to collisions with hydrogen 
atoms, was taken from the precise calculations of Anstee 
& O'Mara (IT995I Barklem & O'Mara (fT997l [1998), and 
Barklem et al. dl998l l. For the other lines, this parameter was 
derived from the fit of the solar atlas (Kurucz et al., 1984 >. 
These values are, for the high excitation lines of Ca II, higher 
than the values obtained from the Unsold or Kurucz approxi- 
mation. However, they agree well with the values obtained by 
Ivanova et al. (2002), who note that the use of the Unsold or 
Kurucz approximation often leads to an underestimation of the 
broadening parameter. 

• We retrieved the spectrum of Procyon from the UVES 
POP (Bagnulo et al.. 12003b . Procyon has a solar-like chemical 
composition, and its surface gravity, derived from Hipparcos 
measurements (Perryman et al., 11997) , is log g=3.96. We 
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Table 2. Parameters of the Ca lines used in the solar spectrum. 
71 = logyvw/Nf{ for a temperature of 10000K. 



line 

(nm) 




Ref 


7i 


Ref 


line 
(nm) 




Ref 


7i 


Ref 


Cal 










Cal 










410.8526 


-0.824 


1 


-6.97 


11 


616.6439 


-1.143 


4 


-7.15 


9 


422.6728 


0.244 


6 


-7.56 


9 


616.9042 


-0.797 


4 


-7.15 


9 


428.3011 


-0.220 


3 


-7.50 


12 


616.9563 


-0.478 


4 


-7.15 


9 


428.9367 


-0.300 


3 


-7.50 


12 


643.9075 


0.390 


4 


-7.57 


9 


430.2528 


0.280 


3 


-7.50 


12 


645.5598 


-1.340 


2 


-7.65 


1 


431.8652 


-0.211 


3 


-7.50 


12 


646.2567 


0.262 


4 


-7.57 


9 


435.5079 


-0.420 


3 


-6.80 


11 


647.1662 


-0.686 


4 


-7.57 


9 


442.5437 


-0.360 


3 


-7.16 


9 


649.3781 


-0.109 


4 


-7.57 


9 


443.4957 


-0.010 


3 


-7.16 


9 


649.9650 


-0.818 


4 


-7.57 


9 


443.5679 


-0.523 


3 


-7.16 


9 


657.2779 


-4.296 


6 


-7.69 


9 


445.4779 


0.260 


3 


-7.16 


9 


671.7681 


-0.523 


4 


-7.14 


9 


445.6616 


-1.660 


3 


-7.16 


9 


714.8150 


0.137 


4 


-7.80 


1 


451.2268 


-1.892 


4 


-7.26 


1 


720.2200 


-0.262 


7 


-7.80 


1 


452.6928 


-0.548 


4 


-7.00 


11 


732.6145 


-0.208 


7 


-7.30 


11 


457.8551 


-0.697 


4 


-7.00 


11 


1034.3820 


-0.409 


5 


-7.12 


9 


468.5268 


-0.880 


3 


-6.40 


1 1 












518.8844 


-0.115 


4 


-7.00 


11 


Call 










526.0387 


-1.719 


4 


-7.51 


9 


393.3663 


0.104 


6 


-7.76 


9 


526.1704 


-0.579 


4 


-7.51 


9 


396.8469 


-0.201 


6 


-7.76 


9 


526.5556 


-0.114 


4 


-7.51 


9 


500.1479 


-0.506 


10 


-7.34 


1 


534.9465 


-0.310 


4 


-7.44 


12 


645.6875 


0.412 


10 


-6.61 


11 


551.2980 


-0.464 


7 


-7.00 


1 1 


820.1722 


0.300 


2 


-7.00 


1 1 


558.1965 


-0.555 


4 


-7.54 


9 


824.8796 


0.570 


2 


-7.00 


11 


558.8749 


0.358 


4 


-7.54 


9 


825.4730 


-0.400 


2 


-7.00 


1 1 


559.0114 


-0.571 


4 


-7.54 


9 


849.8023 


-1.416 


8 


-7.68 


9 


559.4462 


0.097 


4 


-7.54 


9 


854.2091 


-0.463 


8 


-7.68 


9 


585.7451 


0.240 


4 


-7.12 


11 


866.2141 


-0.723 


8 


-7.68 


9 


586.7562 


-1.570 


7 


-7.00 


11 


891.2068 


0.631 


10 


-7.21 


1 1 


610.2723 


-0.770 


4 


-7.19 


9 


892.7356 


0.808 


10 


-7.21 


11 


612.2217 


-0.319 


4 


-7.19 


9 


985.4759 


-0.228 


10 


-7.00 


11 


615.6023 


-2.497 


5 


-7.15 


9 


989.0628 


1.270 


10 


-7.01 


11 


616.1297 


-1.268 


4 


-7.15 


9 


993.1374 


0.051 


2 


-7.00 


11 


616.2173 


-0.090 


3 


-7.19 


9 


1183.8997 


0.290 


2 


-7.36 


11 


616.3755 


-1.286 


7 


-7.15 


9 


1194.9745 


-0.010 


2 


-7.36 


11 


References 
1: VALD 






2:Wiese et al. (1969) 







3: Wiese & Martin (1980) 4: Smith & Raggett (1981) 

5: Kurucz & Bell (1995) 6: Morton (1991) 

7: Smith et al. (1988) 8: Theodosiou (1989) 

9: Anstee & O'Mara (1995), Barklem & O'Mara (1997,1998) 

10: Opacity project 11: fit of the Solar atlas 

12: Unsold formula + fit of the Solar atlas 



adopted a microturbulence =1.8 kms . For this star, 
Mashonkina et al. d2007b adopting T eff = 6510 K (Mashonkina 
et al., 120031) or T eff = 6590 K (Kom et al., I2U031) found a 
subsolar abundance of Ca and a difference of about 0.2 dex 
between the NLTE calcium abundance derived from the 
Cal and Call lines. 

With T e ff=6510 K and our model atom, we derive for Procyon 
a near solar calcium abundance and a much better agreement 
between Ca I and Ca II: log e(Ca) = 6.25 +0.04 from Ca I lines 



and 6.27 +0.06 from Call lines (Fig.[D- 

• We also tested our model on two classical metal-poor 
stars: HD 122563 (a star from our sample of giants) and 
HD 140283. The spectrum of HD 140283 was retrieved from 
the ESO POP (Bagnulo et al., [20031) . For HD 122563 we 
adopted the parameters given in Table [T] and for HD 140283 
the parameters given by Hosford et al. (2009 i: T e tf=5750 K, 
logg=3.4, ft=l-5 kms -1 . The agreement is good. From the 
Cal lines we obtain loge(Ca) = 4.12 + 0.04 and from the 
Call lines log e(Ca) = 4.08 + 0.05. In Fig.|2]we show the fit be- 
tween the observed spectrum of HD 140283 and the synthetic 
profiles computed with log e(Ca) = 4.12. 




4226 4226.5 4227 4454 4454.5 4455 




8660 S662 S664 



Fig. 3. Profiles of three Ca lines computed for a giant star with 
[Ca/H] * -2.5 with LTE (thin blue line) and NLTE (thick red 
line) hypotheses. The wavelengths are given in A. 
a) The NLTE profile of the Cal resonance line is narrower 
in the wings and deeper in the core, b) In a Cal subordinate 
line, for the same abundance of calcium the equivalent width 
computed under the NLTE hypothesis is slightly smaller, c) The 
NLTE correction is important for the strong IR Ca II line (note 
that the scale in wavelength is different for this line), but the 
wings are not affected and a reliable calcium abundance can be 
deduced from these wings via LTE analysis. 



4. Calcium abundance in the EMP stars sample 

4.1. Comparison between the different systems 

In Figj3] we show the influence of the NLTE effects on the 
profile of some typical calcium lines measurable in extremely 
metal-poor stars. Evidently, in particular the wings of the lines 
of the Call infrared triplet are not sensitive to NLTE effects 
until the wings are strong enough. 

In Fig. |4] the correction NLTE-LTE is given as a funtion 
of the temperature, the gravity, and the calcium abundance for 
the Ca I resonance line and a typical subordinate line of Ca I. 
The NLTE corrections for different temperatures, gravities, and 
abundances for all the Ca lines used in this analysis (Table |2]l 
are available in electronic forrrfl The influence of the complex 
NLTE effects as a function of the stellar parameters and of the 

2 http://cdsweb.u-strasbg.fr/cgi-bin/qcat ?J/A+ A/ 
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Fig. 4. NLTE corrections for the Ca I resonance line at 422.67 nm and for a typical subordinate line of Ca I at 445.48 nm. 



characteristics of the lines has been discussed by Mashonkina 
et al. d2007T > and Merle et al. d20TTb . 

In Table Q] we give for each star the mean abundance de- 
rived from the subordinate lines of Ca I (col. 6), the number 
of lines used for the analysis (col. 7), and the internal error 
(col. 8). Column 9 lists the abundance deduced from the res- 
onance line of Cal, and column 10 the abundance computed 
from the infra-red Call line at 866.21nm (only this line of 
the red Call triplet is in the wavelength range of our spectra). 
Columns 11 and 12 list [Ca/H] and [Ca/Fe] deduced from the 
subordinate lines of Ca I. For an easier comparison to the pre- 
vious "First Stars" papers, the solar abundance of calcium was 
taken from Grevesse & Sauval (2000) as in Cayrel et al. (2004 1: 
loge(Ca) © = 6.36. The oscillator strengths of the calcium lines 
measured in our metal-poor stars were updated (to be the same 
as in section [3~2l . As a consequence, for some lines (Table |2]i 
they are slightly different from the log gf values used in Cayrel 
et al. (120041 and Bonifacio et al. (120091 . 

The abundance of Ca was deduced from the equivalent 
widths for the subordinate lines of Cal. But for the reso- 
nance line, which is often strong and slightly blended with Fe I 
and CH lines, a fit of the profile was made. 
The Ca abundance was generally deduced from a fit of the 
wings (insensitive to NLTE effects) for the red Ca II lines. 
However, in the turnoff stars and in the most metal-poor gi- 
ants, the wings almost disappear (when the equivalent widths 
are less than 400 mA) and in this case, equivalent widths were 
used. 

In Fig. [5] we compare the abundances of calcium derived 
from these different systems. The abundance deduced from the 



line of the infrared triplet of Ca II at 862.21nm is, as a mean, 
0.07 dex lower than the abundance deduced from the subor- 
dinate lines of Cal (Fig. [5^). This small difference is quite 
satisfactory since the abundance deduced from the Ca II lines 
depends on the surface gravity of the model, and this gravity 
can be affected by a systematic error since it has been derived 
from the ionization equilibrium of iron under the LTE hypoth- 
esis (see Cayrel et al., 2004). 

In Fig[5j>, the calcium abundance derived from the subordi- 
nate lines of Ca I is compared to the abundance deduced from 
the Cal resonance line at 422.67nm. It is well known that in 
very metal-poor stars, under the LTE hypothesis, this resonance 
line leads to an underestimation of the calcium abundance: e.g. 
Magain (11988b . Ryan et al. d 19961 1. But it was expected that 
NLTE computations of the lines would remove this effect. For 
one typical metal-poor giant CS 22172-02 we show in Fig. [6] 
the observed spectrum and theoretical spectra computed with 
abundances log e(Ca) = 2.52 (best fit), 2.8 and 3.1. The subor- 
dinate lines lead to a Ca abundance of 3.11; this value is well 
established: seven of the subordinate lines have an equivalent 
width stronger than 9 mA and the scatter in abundance is only 
0.09 dex. 

A similar discrepancy has been observed, after NLTE com- 
putations, in several metal-poor stars by Mashonkina et al. 
(2007) and the authors suggested that the explanation could lie 
with the ID atmospheric models adopted for the computations, 
since the Ca I resonance line is formed over a more extended 
range of atmospheric depths than the subordinate lines. 
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Fig. 5. Difference between [Ca/H] deduced from the subordi- 
nate lines of Ca I and [Ca/H] deduced from a) the infrared 
triplet of Ca II and b) the resonance line of Ca I. The filled 
circles represent the turnoff stars and the open symbols the gi- 
ants. The error bar on [Ca/H], A Ca 866.21 and A Ca 422.67 is 
generally less than 0.1 dex. It is indicated only when it exceeds 
0.1 dex. 

a) The Ca abundance deduced from the line of the infrared 
triplet is, as a mean, 0.07 dex lower than the abundance de- 
rived from the subordinate lines of Ca I. 

b) For [Ca/H] tss -2, the abundance deduced from the Ca I res- 
onance line agrees quite well with the abundances deduced 
from subordinate lines, but a discrepancy appears and increases 
linearly when [Ca/H] decreases. It reaches about 0.4 dex for 
[Ca/H]= -3.5. 
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Fig. 6. Observed profile of the resonance line of Ca I (crosses) 
for one typical metal-poor giant compared to theoretical pro- 
files computed with log e(Ca) = 2.52 (thick red line), 2.8 
and 3.1 (thin blue lines). The wavelengths are in A. The (very 
small) difference between the observed spectrum and the pro- 
file computed with log e(Ca) = 2.52 is shown at the bottom of 
the figure (dots, shifted by 0.1). For this star, the subordinate 
lines lead to log e(Ca) = 3.11 (Table [TJ. 

4.2. Influence of convection - 3D models 
4.2.1 . Shift of the calcium lines 

The atmospheres of cool stars are not static. Velocity and in- 
tensity fluctuations caused by convection are observed in the 



Sun (Rutten et al., 2004) and in metal-poor stars. Ramirez et al. 
(2010) found that in HD 122563 (a star of our sample of giants) 
the cores of the Fe I lines are shifted relative to the mean ra- 
dial velocity of the star, this shift increases with the equivalent 
width of the line. Weaker lines form in deeper layers, where the 
granulation velocities and intensity contrast are higher. We also 
observe this phenomenon for the calcium lines in all stars of our 
sample. The radial velocity of the stars were determined from 
a constant set of iron lines, and relative to this "zero point", the 
radial velocity derived from the Ca I resonance line is about 
0.4 kms~' higher in the giants and 0.2 kms 1 higher in the 
dwarfs. In contrast, the radial velocity derived from the (weak) 
subordinate calcium lines is about 0.4 kms~' lower in the gi- 
ants and about 0.2 kms -1 lower in the dwarfs. (We were able 
to measure precisely the shift of the subordinate lines only on 
the "blue spectra" when these lines were larger than 20 mA.) 
At the resolution of our VLT/UVES spectra with R » 45000, 
the asymmetry of the lines cannot be reliably measured. 

The shift of the Ca I resonance line does not clearly depend 
on [Ca/H], and therefore it does not seem that there is a clear 
correlation between the shift of the Ca lines and the discrep- 
ancy between the abundances deduced from the resonance or 
the subordinate lines. 

4.2.2. Abundance correction 

The largest abundance corrections caused by granulation ef- 
fects occur at low metallicities. This is mainly because the dif- 
ference between the ID and 3D predictions for the mean tem- 
perature of the outer layers of metal-poor stars is very larg e (see 
e.g. Gonzalez Hernandez et al., 1201 Oi l. We tried to investigate 
the change in abundances caused by thermal inhomogeneities 
and differences in formation depth (3D corrections hereafter) 
for the Ca I resonance line at 422.67 nm and for a typical sub- 
ordinate line of Ca I at 445.48 nm. 

-For a representative tumoff star, we used a 3D-C05BOLD 
model (Freytag et al., 120021 [2QTTb from the CIFIST grid (see 
Ludwig et al., 12009b with parameters (T e ff, logg, [Fe/H]): 
6270K/4.0/-3.0). 

-For the giants we used two models (4488 K/2.0/-3.0 and 
5020 K/2.5/-3.0) from the CIFIST grid. They have both a grav- 
ity slightly higher than the ones in our sample of giants, but no 
closer 3D model is available at the moment. 

From these computations we found that for [Fe/H] =-3, 
[Ca/H] =-2.6, the 3D correction for the subordinate lines of 
Cal in turnoff and in giant stars is very small. 
In giants, the 3D correction seems to be negligible also for the 
resonance line. But with the model of turnoff stars, we found 
for the resonance line a 3D correction of -0.44 dex, and there- 
fore the 3D correction increases the discrepancy between the 
subordinate lines and the 422.67 resonance line. 

To date, it is not well understood why the abundance of 
calcium deduced from the resonance line of Ca I is at very low 
metallicity, lower than the abundance deduced from the subor- 
dinate lines. It would be interesting to repeat the 3D compu- 
tations with more metal-poor models (not available today). A 
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fully consistent 3D NLTE model atmosphere and line forma- 
tion scheme is currently beyond our reach. 

5. Results and discussion 




Fig. 7. [Ca/Fe] vs. [Fe/H] in our sample of EMP stars. Symbols 
as in Fig [5] Turnoff and giant stars agrees quite well, [Ca/Fe] 
is constant in the range -4.3 < [Fe/H] < -2.5 and the mean 
value is equal to [Ca/Fe] =0.5, a value slightly higher than the 
LTE value [Ca/Fe]=0.35 (see Bonifacio et al., |2"0"09"l >. 
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Fig. 8. Abundance ratios of O, Na, Mg, Al, S, and K relative to 
Ca in the early Galaxy. The abundances of all these elements 
were computed taking into account the NLTE effect. The black 
dots represent the turnoff stars, the open circles the giants (blue 
for the unmixed giants and red for the mixed giants). 



In this section we adopt the Ca abundance deduced from 
the subordinate lines of Cal (Tab.Q] columns 11 and 12). 

In Fig. [7] we present the new relation between [Ca/Fe] and 
[Fe/H]. The error bar plotted in for the subordinate lines the fig- 
ure is the quadratic sum of the error due to the uncertainty of the 
model and the random error of the Ca abundance derived from 
the subordinate lines. The agreement between the turnoff and 
the giant stars is excellent, the ratio [Ca/Fe] is constant in the 
interval -4.5 < [Fe/H] < -2.5. The slope of the regression line 
is -0.05. The scatter of [Ca/Fe] is a little smaller when NLTE 
effects are taken into account: 0.09 from NLTE computations 
compared to 0. 10 from LTE computation . 

In Fig. [8] we present the behavior of the abundances of O, 
Na, Mg, Al, S and K relative to Ca. The computation of these 
abundances take into account the NLTE effects (Andrievsky et 
al. 120071 120081 [20101 and Spite et al. I20TTI Since the abun- 
dance of oxygen has been determined from the forbidden oxy- 
gen line at 630 nm (Cayrel et al. 2004), it is free of NLTE ef- 
fects. 

In Fig. [8] different symbols are used for mixed and unmixed 
giants. After Spite et al. (2005, 2006) we call "mixed giants", 
those where, ownng to mixing with deep layers, carbon has 
been partially transformed into nitrogen ([C/N] < -0.6), part 
of 12 C has been transformed into 13 C ( 12 C/ 13 C < 10), and 
lithium is not detected (lithium has been severely depleted by 
this mixing). In the HR diagram, these "mixed giants" are lo- 
cated above the "bump". It has been found (Andrievsky et al., 
2007) that some mixed giants are enriched in sodium, this is 
visible in Fig.[8j). This Na-enhancement reflects a deep internal 
mixing (or an AGB status, or a contamination by AGB stars), 
but it does not reflect an anomaly of the chemical composition 
of the cloud that formed the star. Therefore the mixed stars c an- 
not be used to determine the ratio [Na/Ca] in the early galactic 
matter. 



Table 3. Mean value of the ratios [X/Fe], [X/Ca] and [X/Mg] 
in the interval -3.6 < [Fe/H] < -2.5 and scatter around the 
mean 
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Since the abundance ratios of O, Na, Al, S, and K rela- 
tive to Ca are fairly flat vs. [Fe/H] in the central part of the 
diagram, we can define a mean value in this central interval, 
say -3.6 < [Fe/H] < -2.5 as has been done with Mg in 
Andrievsky et al. (BUTOl l. The mean values of [X/Fe], [X/Mg] 
and [X/Ca] are given in Table|3]with the corresponding scatter. 
However, an anti-correlation seems to exist between [S/Ca] and 
[Fe/H] and also [S/Mg] and [Fe/H]: the Kendall r coefficient is 
98.9% for [S/Ca] and 99.1% for [S/Mg]. In these computations 
the weight of BD +17°3248 ([Fe/H]=-2.2) is important. This 
star is quite peculiar: according to For & Sneden (12010) it is a 
red horizontal branch star strongly enriched in heavy elements 
(Cowan et al., 2002). However, if this star is removed from the 
computations, the Kendall r coefficient remains high: 97.6% 
for [S/Ca] and 97.9% for [S/Mg] . 

One turnoff star, BS 16023-046, and one unmixed giant 
CS 22956-50, have abnormally strong and broad sodium D 
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lines. Their radial velocity (-7.5 kms _I and -0.1 kms -1 ) is low, 
and the stellar lines are very probably blended with interstellar 
lines. These stars have not been plotted on Fig.[HJ) and were not 
taken into account in computing the mean value and the scatter 
of [Na/Fe], [Na/Mg] and [Na/Ca] in Tabled 

It is interesting to take advantage of the high quality of the 
data to compare in Table [3] the scatter around the mean value 
when Fe, Mg, or Ca are taken as reference elements. The scatter 
of [O/Fe], [O/Mg] or [O/Ca] is fairly large: about 0.2 dex but 
this can be explained by the scatter of the oxygen abundance 
due to the reduced number of oxygen lines and the difficulty 
of the measurements. For the other ratios the scatter lies be- 
tween 0.09 and 0.15 dex and we consider that this difference 
of 0.06 dex is significant. The scatter is always significantly 
smaller when Ca is taken as the reference element, with two 
exceptions: 

• Oxygen: but in this case we have seen that the scatter is domi- 
nated by the error on the measurement of the very weak oxygen 
line. 

• Aluminum: Al is better correlated with Mg than with Ca. But 
this correlation (opposite to the well known anti-correlation 
found in globular cluster stars) has at least one exception: 
CS 29516-024 is rather Al-poor but it is (relatively) Ca-rich 
and Mg-rich, and it does not seem possible to explain these 
differences by errors of measurements. A similar correlation 
between [Al/Fe] and [Mg/Fe] has been also suggested by Suda 
et al. d201 11 1 even in a large but very inhomogeneous sample of 
metal-poor stars. 

6. Concluding remarks 

We determined the calcium abundance in a homogeneous sam- 
ple of of 53 metal-poor stars (31 of them with [Fe/H] < -2.9) 
taking into account departures from LTE. We have shown that 
the trend of the ratio [Ca/Fe] vs. [Fe/H], below [Fe/H] =-2.7 
is almost flat and derived a new mean value of [Ca/Fe] in the 
early Galaxy: [Ca/Fe] = 0.5 + 0.09. 

Generally speaking, our NLTE calculations agree quite well 
with those of Mashonkina et al. (2007). However, Mashonkina 
et al. had found different Ca abundances from the Cal and 
Call lines in Procyon. With our new model atom of cal- 
cium, we are able to reconcile the calcium abundance deduced 
from neutral and ionized calcium lines in Procyon. We derived 
loge(Ca) = 6.25 +0.04 from Cal lines and 6.27 +0.06 from 
Ca II lines (6.33 and 6.40 under the LTE hypothesis). Moreover, 
the calcium abundance is found to be solar, as expected. 

In metal-poor stars (below [Ca/H]=-2.5), a discrepancy 
clearly appears between the Ca abundances deduced from ei- 
ther the resonance Ca I line, or the Ca I subordinate lines. A 
rough estimation of the effect of convection on the profile of 
these lines does not explain this discrepancy. 

In the stars of our sample (giants and turnoff stars) the 
wavelengths of the calcium lines are shifted by convection as a 
function of the equivalent width of the lines as has been found 
for the iron lines by Ramirez et al. (2010i in HD 122563. 

The scatter around the mean value of [Ca/Fe] is small but 
since the NLTE correction is about the same for all the subor- 
dinate lines used in this analysis, the scatter is almost the same 



as it was in the LTE analyses (Cayrel et al., 2004, Bonifacio et 
al.. l2009l . 

The abundance of the light metals O, Na, Al, S, K is well 
correlated with the calcium abundance. The correlation is gen- 
erally a little better than it is with the magnesium abundance 
(with exception of the tight aluminium/magnesium correla- 
tion). However, it is striking that in Table|3] even when NLTE is 
taken into account, the correlation of the light elements O, Na, 
Al, S, K with Ca or Mg (all supposed to be formed mainly by 
hydrostatic fusion of C, Ne or O), is never better than the corre- 
lation with iron, providing some support to supernovae models 
predicting nucleosynthesis of all these elements predominantly 
in explosive modes (e.g. Nomoto et al. 2006). 
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